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Abstract 

The recently proposed Fully-Renormalized QRPA (FR-QRPA), which fullfils 
the Ikeda sum rule (ISR) exactly, is applied to the two-neutrino double beta 
decay of 76 Ge, 82 Se, 100 Mo, 116 Cd, 128 Te and 130 Xe. The results obtained are 
compared with those of other approaches, standard QRPA and self-consistent 
QRPA (SCQRPA). The similarities and the differences among the methods 
are discussed. The influence of the restoration of the Ikeda sum rule on the 
2^/3/3-decay amplitude is analyzed. 
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I. INTRODUCTION 



Observation of the neutrinoless double beta decay (0z//3/3-decay), violating the total lep- 
ton number by two units, would give unambiguous evidence for new physics beyond the 
Standard Model [1-3]. For instance, at least one of the neutrinos would have to be a Majo- 
rana particle with non-zero mass [4]. The current experimental upper limits on the Ovfifi- 
decay half-life impose stringent constraints, e.g., on the parameters of Grand Unification 
and super-symmetric extensions of the Standard Model. 

Rates of the 2z//3/3-decay, which is a second order process allowed within the Standard 
Model, can be calculated within the same nuclear structure models. Thus, the results of the 
nuclear structure calculations can be directly compared with the corresponding experimental 
data available for a number of nuclei [6]. Such a comparison provides a very useful test of 
the models. 

The Quasiparticle Random Phase Approximation (QRPA) [5] has been successfully ex- 
ploited in nuclear physics to describe properties of the excited states of open-shell nuclei 
and to calculate intensities of various nuclear reactions, including the double beta decay 
(see reviews [1]). 

It was shown that the experimental data on the 2z//5/9-decay rates can be reproduced in 
QRPA calculations with a sufficiently large strength of the particle-particle interaction [7] . 
But the proximity of the value to the point of the QRPA collapse questions the reliability 
of the results. It is known that QRPA collapse occurs due to the use of the quasi-boson 
approximation (QBA) which violates the Pauli exclusion principle (PEP) and generates too 
many ground state correlations. 

Renormalized QRPA (RQRPA) was formulated in Refs. [8] to restore PEP in an ap- 
proximate way. The main goal of the method is to use a self-consistent iteration of the 
QRPA equation with taking into account quasiparticle occupation numbers in the QRPA 
ground state. That leads to a modification of the commutation relations for bifermionic 
operators as compared to the ordinary quasiboson approximation (QBA). At the same time 
so-called scattering terms (describing transitions of the quasiparticles) are neglected in the 
Hamiltonian and in the phonon operators. The RQRPA does not collapse for physical values 
of the particle-particle interaction strength and has been extensively used to calculate the 
intensities of the double beta decay [1,10,11]. It has been also shown that the RQRPA pro- 
vides better agreement with the exact solution of the many-body problem within schematic 
models, even beyond the critical point of the standard QRPA (see, e.g. [9] and references 
therein) . 

The self-consistent RQRPA (SCQRPA) is a more complex version of RQRPA to describe 
the strongly correlated Fermi systems. Within this method one goes a step further beyond 
the RQRPA. In the SCQRPA at the same time the quasiparticle mean field is changed by 
minimizing the energy and fixing the number of particles in the correlated ground state 
of RQRPA instead of the uncorrelated one of BCS as is done in the other versions of the 
RQRPA. In this way SCQRPA partially overcomes the inconsistency between RQRPA and 
the BCS approach and is closer to a fully variational theory. 

Nevertheless, the main drawback of the modern versions of RQRPA and SCQRPA is 
the violation of the model-independent Ikeda sum rule (ISR) [11-13]. A modification of the 
phonon operator by including scattering terms is needed in order to restore the ISR within 
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RQRPA. The fully-Renormalized QRPA (FR-QRPA) was formulated in Ref. [16] for even- 
even nuclei in such a way that it complies with restrictions imposed by the commutativity 
of the phonon creation operator with the total particle number operator. It was shown 
analytically that the Ikeda sum rule is fulfilled within the FR-QRPA [16]. Also FR-QRPA 
is free from the spurious low-energy solutions which would be generated by the scattering 
terms considered as additional degrees of freedom as suggested in [17]. 

The aim of the paper is twofold. First, we would like to describe the FR-QRPA equations 
in more details (as compared with the original paper [16]) for a simple case of a Hamiltonian 
with the separable residual interaction in both particle-hole and particle-particle channels. 
Second, the first numerical application of FR-QRPA is given to calculate 2u(3(3-decay in- 
tensities and relevant quantities. So far, the full convergence of the FR-QRPA solution has 
been obtained only for a rather small model space. Nevertheless a comparison of the results 
obtained within FR-QRPA and SCQRPA can be provided. 



II. BASIC RELATIONSHIPS OF THE FULLY-RENORMALIZED QRPA 



Within RPA an excited nuclear state, with angular momentum J and projection M, is 
created by applying the phonon operator Q JM to the vacuum state |0^ Pj4 ) of the initial, 
even-even, nucleus: 



| JM) = Q\ M K PA ) with QjmKpa) = 0. 



(1) 



As was shown in Ref. [16], the most appropriate way is to write down the phonon 
structure in terms of the particle creation and annihilation operators. That allows to fulfill 
the important principle of the commutativity of Q \ M with the total particle number operator 
A = N + Z. The phonon operator has the following structure: 



QjM 



[ x { P n,j)C\pn, JM)- y<j, n ,j)C(pn, JM) 

pn 



(2) 



with C*(pn,JM) = [c],c n ] JM and C{pn,JM) = {-) J ~ M C{pn,J - M), where c+ mr (c rmT ) 
denotes the particle creation (annihilation) operator for protons and neutrons (r = p,n). 

Going into the quasiparticle representation, the quasiparticle creation and annihilation 
operators a+ mr and a TrrtT , (r = p, n) can be defined by the Bogolyubov transformation 



that leads to the following expression for the phonon operator Q\ M 
Q\m = E \x {pnJ) A^pn, JM) - Y (pnJ) A (pn, JM) 



pn 



(3) 



(4) 



A 1 " = A^ + (u n v n Bi - UpVpB 



(*>n - V p) 



A\pn, JM) = 



44 



JM 



; B^(pn,JM) = 



dp(l n 



JM 
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where X = u p v n x — v p u n y, Y = u p v n y — v p u n x. The bifermionic operators A>, A now being 
the basic building blocks of the FR-QRPA automatically contain the quasiparticle scattering 
terms which, however, are not associated with any additional degree of freedom. That means 
that there are no spurious low-lying solutions in the present theoretical scheme which would 
be generated by the scattering terms considered as independent constituents of the phonon 
operator (as proposed in [17]). 

From this point we can follow the usual way to formulate the RQRPA [8], substituting 
A by A everywhere. The forward- and backward-going free variational amplitudes X and Y 
satisfy the equation: 



A B 
B A 



£ 



U 
-U 



(5) 



where m marks different roots of the QRPA equations for a given J n , 



A=(0^ PA 
B = -Kpa 



A. 



A 



H, A f 
H, A 
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RPA1 1 



|0rpa)> 



and the renormalization matrix U pn is 



Upn = K PA \ \A(pn, JM),A\p'n', JM)] \0 + RPA ) = 6 pp/ 6 nn ,V, 



(6) 



(7) 



We use a rather simple, but realistic, Hamiltonian H consisting of the quasiparticle mean 
field H and the residual separable particle-hole (ph) and particle-particle (pp) interactions: 



H 



H + H* + H» 



inti 



Y, E T a\a T , 



r=p,n 



Hint — XE( 1) M (PimPi-M + Pi-M@1m)i 
M 

Hint = ~ K Yl{ — {PlM^l-M + A-mAm)i 
M 



(8) 
(9) 

(10) 
(11) 



with f3 1M = —J 1 T,(p\\cr\\n) 



pn 



if s. 
-< p W, 



1M' 



P 1M = J' 1 E(p\W\\n) 



pn 



^p^n 



1M 



and J = 1. 



Taking into account the exact (fermionic) expressions for the commutators in (6), (7), 
one gets the following expressions for the FR-QRPA matrices A and B: 



A= (E p + E n )V pn - 2(E P - E n )(uivt + <<)ft ; 



^pn 



Spp' $nn' 



+2x{u p v n u p iv n i + v p u n v p iu n i)V pn V p i n , 
-2K{u p u n u p iu n it> pn V p , n i + v p v n v p iv n iV pn V p i n i) , 



(12) 



13 Qi^Ep E n ^UpVpU n v n T^p n 6ppf6 nn ' 

+2x{u p v n v p m n i + v p u n v p iu n i)V pn V p i n , 

+2K(u p u n v pl v nl V pn t> plnl + v p v n u p <u n it) pn t) p i n <) . (13) 

The renormalization matrices T>,T>,T> entering (7), (12), (13) can be represented in terms 
of the relative quasiparticle occupation numbers Af p for the level p in the RQRPA vacuum: 
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V pn = 1 - M n - M p + (l - v 2 p - vl) K pn) ( 14) 

Vpn = 1 - N n - U p - (u 2 p + 1l pn , 

V pn = l-M n -M p + (y 2 p + v 2 n ) Tl pn , 
with lZ pn = ^J^sT ■ In turn, the quasiparticle occupation numbers 

K = Jr 2 ^RPA\J2 a lm T arm T \0^ pA ); T=p,U. (15) 

can be expressed in terms of the backgoing amplitudes Y of the RQRPA solution (5) [8]. 
In the calculation we shall use the aproximate expression for J\f p and J\f n : 



n \J,m 



pn 



K « V E E( 2 ^ + VKJ 2 v pn (i6) 

V \J,m ) 

where j = y/2j + 1. In the present paper we consider only J 71 " = 1 + contribution to the 
sums in (16). Along with the modified SCQRPA and FR-QRPA equations for the chemical 
potential: 

KpaWK PA ) = Ein On + K - vl)N n ) = N, 
n 

Kpa\ZKpa) = Ei 2 (4 + ("J " OK) = Z, (17) 
p 

a rather complicated set of equations (5)-(17) has to be solved. 

It is noteworthy that the renormalization matrices (14) become the same, T> pn = T> pn = 
T> pn , in the limit TZ pn = and coinciding with the renormalization matrix of the usual 
RQRPA (see, e.g., [11]). Thus, one can argue that the standard versions of RQRPA neglect 
effectively the differences between the quasiparticle occupation numbers whereas SCQRPA 
and FR-QRPA take the differences into account. 

From now on we follow the usual way of solving RQRPA equations [8]. It is useful to 
introduce the notation: 

X = U 1/2 X, Y = U 1/2 Y, (18) 

A = U- 1 ' 2 AU- 1 / 2 , B=U- l ' 2 BU- 1 ' 2 . (19) 
Then the amplitudes X and Y satisfy the equation of usual QRPA: 

{& a) (? m ) =Sm (o -i) ' (20) 

Solving the FR-QRPA equations, one gets the fully renormalized amplitudes X, Y with the 
usual normalization and closure relations: 
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Evm yfc vm \rk _ r 

A (pn,J) A (pn,J) - Y (pn,J) Y (pn,J) ~ °km, 

pn 

^(pnJ^ip^uJ) 1 (jmj) 1 (p^uJ) - Vpp^nm, 

m 

A (pf»,J) r (p 1 ni,J) ~ r (pn,J)^(p ini ,J) ~ U - 

m 

It was shown analytically that the Ikeda sum rule is fulfilled within the FR-QRPA [16], 
in contrast to the earlier versions of the RQRPA [8] . The Ikeda sum rule states that the dif- 
ference between the total Gamow- Teller strengths and in the (3~ and (3 + channels, 
respectively, is 3(N - Z) [14]: 

ISR = SM -S™ = 3(N-Z), (22) 



Mm 



S<->=£ (l + M,m\fc M K PA ) , 5(+)= E (l + M,m'|/3+ M |0+ 



Mm' 



(23) 



With the use of the closure conditions (21), the expressions for T> pn (14) and the chemical 
potentials (17), one can show [16] that 

ISR = £ \(p\\qj\\n)\ 2 (v 2 n - v 2 p )V pn = 3(N - Z). (24) 

pn 

The inverse half-life of the 2z//3/3-decay can be expressed as a product of an accurately 
known phase-space factor G 2v and the second order Gamow- Teller transition matrix element 

KMs. - 0+J]- 1 = G 2u (g A r \M 2 J T \ 2 . (25) 

The contribution from the two successive Fermi transitions is safely neglected as they arise 
from isospin mixing effect [2]. The double Gamow- Teller matrix element M"q T for ground 
state to ground state 2z//3 / 3-decay transition acquires the form 

M2 „ = v (Q/ + ll P- II i^Xi^li^Xi^ II P- II o, + ) 

The sum extends over all 1 + states of the intermediate nucleus. The index i(f) indicates that 
the quasiparticles and the excited states of the nucleus are defined with respect to the initial 
(final) nuclear ground state |0+) (|0j)). The overlap is necessary since these intermediate 
states are not orthogonal to each other. The two sets of intermediate nuclear states generated 
from the initial and final ground states are not identical within the considered approximation 
scheme. Therefore the overlap factor of these states is introduced in the theory as follows: 

= E[^n(l + m i )X p „(l + m / ) - Y pn {l+ mi )Y pn {l + m f )]. (27) 

pn 
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III. CALCULATION RESULTS 



In this section we present the 2z//3/3-decay results obtained within the FR-QRPA for 
76 Ge, 82 Se, 100 Mo, 116 Cd, 128 Te and 130 Xe, in comparison with the QRPA and SCQRPA ones. 
Rather small model bases listed in the Table I are used in order to get full convergence in the 
FR-QRPA method. The levels are in a vicinity of the Fermi levels and spin-orbit partners are 
always taken into account. FR-QRPA method is rather sensitive to the differences between 
occupation probabilities for protons and neutrons entering the denominator in the expression 
of the bifermionic operators A* , A (4) and in the expresion for TZ pn factor of renormalization 
matrices (14). 

For levels far from the Fermi one, the values for occupation probabilities for protons 
and neutrons become almost equal. Because of the denominator which appears in the 
expression of bifermionic operators (4), that causes numerical problems, in particular the 
method doesn't converge for large enough values of particle-particle strength. Therefore, the 
bases are fixed in order to get convergence for a larger interval of particle-particle strength, 
in particular up to the point of the collapse of the 2z/ / 3 / 3-decay matrix elements. 

The single particle energies are obtained by using a Coulomb-corrected Woods-Saxon 
potential with Bertsch parametrization. The proton and neutron pairing gaps are determined 
phenomenologically to reproduce the odd-even mass differences through a symmetric five- 
term formula [18]. Then the equations for the chemical potentials (17) are solved for proton 
and neutron subsystems. The pairing gaps entering the BCS equations are given in the 
Table I. 

The calculation of the QRPA energies and wave functions requires the knowledge of 
the particle-hole x an d particle-particle k strengths of the residual interaction. The value 
of particle-hole strength x parameter for each nucleus is fixed in order to reproduce the 
experimental position of the Gamow- Teller giant resonance in odd-odd intermediate nucleus 
as obtained from the (p,n) reactions [19], [20], [21]. Those values are also given in the Table 
I. The particle-particle strength k is considered as a free parameter. 

The numerical results are shown for two groups of nuclei, the nuclei with A < 100 and 
A > 100 respectively. The calculations are done within QRPA, SCQRPA and FR-QRPA in 
order to show the better stability of the latter method. 

The main drawback of the QRPA is the overestimation of the ground state correlations 
leading to the collapse of the QRPA ground state, near a certain critical interaction strength. 
Around this point the backward-going RPA amplitudes Y pn of the first 1 + states become 
overrated, and too many correlations in the ground state are generated with increasing 
strength of the particle-particle interaction. This phenomenon, as a result of the quasiboson 
approximation used, leads to QRPA collapse and implies an ambiguous determination of 
the (3 and 2z//9/9-decay matrix elements. 

In Fig.l and Fig.2 the dependence of the energy of the first excited Gamow- Teller state in 
daughter nuclei is plotted versus the k, parameter. Hereafter, the dashed line corresponds to 
the QRPA case, the dotted line represents the SCQRPA case and the solid line describes FR- 
QRPA case. For all studied nuclei the collapse of the first excited state is shifted to higher 
values of k for each method and the stability increases in the FR-QRPA case. In Fig. 3 and 
Fig. 4 the 2z/ / 3 / 3-decay matrix elements as a function of the particle-particle strength k are 
shown. The calculations are done for all nuclei within the three metods. The horizontal 
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dashed line indicates the experimental values taken from [3]. 

For all nuclei there is a similar behaviour in the sense that QRPA and SCQRPA collapse 
a bit earlier than the FR-QRPA does. Although the chosen bases are rather small, the new 
effects we intend to emphasise as the differences between QRPA extensions are evident. The 
FR-QRPA method offers considerably less sensitive dependence of M"q T on k and shifts the 
collapse to larger values of particle- particle strength. 

According to the definition (22) S~ {S + ) is the total summed Gamow- Teller (3~{f3 + ) 
transition strength from the ground state of an even-even nucleus. In Fig. 5 and Fig. 6 we 
plot the relative (3~ strength, S~/3(N — Z) for the mother nucleus (left side) and relative (3 + 
strength, S + /3(N - Z) for the daugther (right side), for A < 100 and A > 100 respectively 
as a function of particle-particle interaction parameter in order to show the magnitude and 
the nature of violation of ISR. 

Now we would like to discuss the conservation of the Ikeda sum rule ISR = S_ — S + = 
3(N — Z) in the FR-QRPA framework and to compare with the previous calculations for 
QRPA ans SCQRPA. We didn't include the calculations for RQRPA because SCQRPA goes 
beyond and brings more improvements than RQRPA, especially for Ikeda sum rule. 

Finally we combine the data of previous plots and show the ratio of SCQRPA and FR- 
QRPA sum ISR/3(N — Z) as a function of k. In the QRPA the Ikeda sum rule is exactely 
conserved as long as all spin-orbit partners of the single-particle orbitals are included. In 
the other extended versions of QRPA the sum rule is violated with a degree of deviation 
lying between f 7% (RQRPA) and 3% (SCQRPA) [12]. In our study, following the analytical 
calculation of [16], we have shown numerically that Ikeda sum rule is exactly fullfiled within 
FR-QRPA formalism. 

IV. CONCLUSIONS 

In summary, the first calculation of the 2z/ / 9 / 9-decay matrix elements within the recently 
proposed Fully-Renormalized QRPA, which fulfills Ikeda sum rule exactly, are presented. 
The considered nuclear model includes the separable Gamow- Teller residual interaction. 
The subject of interest is the effect of the restoration of the Ikeda sum rule on the 2z//3/3- 
decay observable for A = 76,82,100,116,128,130 systems. Within the present work we 
arrived to the folowing important conclusions: 

i) The SCQRPA violates the Ikeda sum rule. This phenomenon has been indicated in 
the previous studies [12], but the degree of violation we obtained is less than in the other 
calculations because we did not include all multipolarities. 

ii) In the limit when the difference between proton and neutron quasiparticle occupation 
numbers is neglected the FR-QRPA coincides with SCQRPA. 

iii) From a comparison of FR-QRPA with SCQRPA results we conclude that the effect 
of the restoration of the Ikeda sum rule is important in the range of large value of particle- 
particle strength beyond the point of collapse of the standard QRPA. 

It is worth to mention that the FR-QRPA approach is sensitive to the precise evaluation 
of the proton and neutron quasiparticle occupation numbers. Due to the limitation of the 
approximate expression given in (16) (motivated by a similarity to the SQRPA approach) 
the convergence of the FR-QRPA is achieved only for relatively small model space. However, 
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even for such a model space the differences among the standard QRPA, SCQRPA, FR-QRPA 
approach are evident for k close to the point the standard QRPA breaks up. There is a hope 
that for a proper ansatz of the RPA ground state the FR-QRPA approach can work also for 
a large model space. This is the subject of our further study. 
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TABLES 



TABLE I. The proton and neutron pairing gaps determined phenomenologically to reproduce 
the odd-even mass difference and the particle-hole strength \ chosen to reproduce the experimental 
position of Gamow- Teller resonance. The single particle basis for all nuclei under consideration is 
also shown. 





76 Ge 76 Se 


82 Se 


82 Kr 


100 Mo 


100 Ru 


Basis 


lp, Of, Og 


lp, Of, Og 


lp, 2s, Id, Og 


A p [MeV] 


1.561 


1.751 


1.401 


1.734 


1.612 


1.548 


A n [MeV] 


1.535 


1.710 


1.544 


1.644 


1.358 


1.296 


X [MeV] 


0.21 


0.18 


0.17 




116 Cd 


116 Sn 


128 Te 


128 Xe 


iso Te 


130 Xe 


Basis 


lp, 2s, Id, Og, Oh 


2s, Id, Og, Oh 
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FIG. 1. The energy of the first excited Gamow- Teller state as function of particle-particle 
interaction strength k for the daughter nuclei with A < 100 
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FIG. 2. The energy of the first excited Gamow- Teller state as function of particle-particle 
interaction strength k for the daughter nuclei with A > 100 
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FIG. 3. 2^/3/3-decay matrix elements as a function of particle-particle strength k for A < 100 
nuclei 
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FIG. 4. 2z//3/3-decay matrix elements as a function of particle-particle strength k for nuclei with 
A > 100 
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FIG. 5. The dependence of the total f3 and (3 + strength on the k parameter for mother and 
daughter nuclei respectively. The total strength is normalized to 3(N-Z) 
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FIG. 6. The dependence of the total f3 and (3 + strength on the k parameter for mother and 
daughter nuclei respectively. The total strength is normalized to 3(N-Z) 
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FIG. 7. Ikeda sum rule as function of particle-particle interaction strength k nuclei with 
A < 100. The dashed line coresponds to SCQRPA and the continuous one coresponds to FRQRPA. 



19 



1.01 



N 

I 

Z 

X 1 



0.99 
1.01 



i i i i I i i i i I i i i i I i i i i 



1.01 



N 

- i 

^ 1 

J 22 



1 1 1 1 1 1 1 1 1 1 1 

- 116 Sn ' 

- 


II II 

- 


N 

\ 

\ 




\ 

\ 

\ 

\ 

\ 

I I I I I I I I I h I I I I 


I I I I 




0.99 



0.99 



0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 
K[MeV] K[MeV] 



FIG. 8. Ikeda sum rule as function of particle-particle interaction strength k nuclei with 
A < 100. The dashed line coresponds to SCQRPA and the continuous one coresponds to FRQRPA 
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